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ABSTRACT 

We study the incidence rate of damped Lya systems associated with the host galaxies of gamma- 
ray bursts (GRB-host-DLAs) as functions of neutral hydrogen column density (JVhi) and projected 
star formation rate (SFR) using cosmological SPH simulations. Assuming that the occurrence of 
GRBs is correlated with the local SFR, we find that the median A^hi of GRB-host-DLAs progressively 
shifts to lower AThi values with increasing redshift, and the incidence rate of GRB-host-DLAs with 
log Am > 21.0 decreases rapidly at z > 6. Our results suggest that the likelihood of observing the 
signature of IGM attenuation in GRB afterglows increases towards higher redshift, because it will not 
be blocked by the red damping wing of DLAs in the GRB host galaxies. This enhances the prospects 
of using high-redshift GRBs to probe the reionization history of the Universe. The overall incidence 
rate of GRB-host-DLAs decreases monotonically with increasing redshift, whereas that of QSO-DLAs 
increases up to z = 6. A measurement of the difference between the two incidence rates would enable 
an estimation of the value of ??grb , which is the mass fraction of stars that become GRBs for a given 
amount of star formation. Our predictions can be tested by upcoming high-z GRB missions, including 
JANUS (Joint Astrophysics Nascent Universe Scout) and SVOM (Space multi-band Variable Object 
Monitor). 

Subject headings: cosmology: theory — stars: formation — galaxies: evolution - galaxies: formation 
- methods: numerical 



1. INTRODUCTION 

A number of authors have proposed using GRBs 
to probe the history of cosmic star for mation 
and t he reionization of th e Universe (e.g., iTotanil 
19971: iMiralda-Escudd Il998t lLamb fc Reichartl 120001 : 
Barkana fc Loebl 120041). neither of which is well- 



understood (see, e.g.. lHolder et a l. 2003; Nagamin e et al.l 
2006). To date, observations of high-redshift (hereafter 
high-z) quasars and galaxies ha ve been able to c onstrain 
reionization only up to z ~ 7 (|Fan et al.l I2006T ) . How- 
ever, if GR Bs are associated with the death s of massive 
stars (e.g., IWooslevi [l99l iPaczvnskil [l99ll . then theo- 
retical studies imply that GRBs may be detectable out 
to z ~ 10 — 2 through their prompt 7-ray emission 
and a fterglows (lLamb fc Reichartl 120001: ICiardi fc Loebl 
120001 IGou et al l 12004 llnoue et al.l 120071 ). This raises 
the possibility of using GRBs to investigate the reion- 
ization history of the Universe, and the Swift satellite 
has indeed detected high-z GRBs with bright afterglows 
(ICusumano et al.ll200a lHaislip et all l2006t iKawai et all 
2006). 

However, GRB lines-of-sight (LOSs) tend to probe 
more the inner parts of galaxies than the random 
QSO LOSs do, and are therefore often associated 
with neutral hydrogen (Hi) absorption (jProchaska et al.l 
120071) . Indeed, analyses of the afterglow spectra re- 
yeal the presence of DLAs in the red damping wing 
(|Vreeswiik et al.l 12004 iBerger et aT1l2006l : IWatson et all 
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120061: iRuiz-Velasco et al.l 120071 ). and the DLAs may 
hide the absorptio n signatures of the neutral IGM 
(jTotani et all [2006) . If such cases dominate the high- 
z GRB afterglow spectra, then it may be difficult to use 
GRBs to probe the detai l ed re ionization history of the 
Universe (|McQuinn et a l. 2008). Thus, it is important 
to understand the redshift evolution of the incidence rate 
of GRB-host-DLAs at z > 6 as a function of A^hi- 

In this Letter, we use cosmological SPH simulations 
based on the concordance A cold dark matter (CDM) 
model to study the Am distribution and the incidence 
rate of GRB-host-DLAs as a function of redshift between 
z = 1 — 10. The DLAs associated with quasar LOSs are 
often referred to as QSO-DLAs. Since quasars serve as 
randomly distributed background beacons in the Uni- 
verse, QSO-DLAs can be more broadly interpreted as 
all the Hi gas clouds that satisfy the DLA criterion 
(AThi > 2 x 10 20 cm -2 ), regardless of whether or not they 
have been intersected by quasar LOSs. We adopt the lat- 
ter broad interpretation of QSO-DLAs in this paper, and 
by this definition GRB-host-DLAs are a subset of QSO- 
DLAs. 

2. SIMULATIONS 

Our simulations were performed with the smoothed 
particle hydr odynamics (SPH) code GADGET-2 
(Springel 2005), which includes radiative cooling by hy- 
drogen and helium, heating by a uni form UV background 
(|Katz et alj|1996t iDave et al.lll999f ). star formation and 
supernova feedback based on a sub-particle multiphase 
ISM model, and a phenomenological de scription of 
galactic winds (|Springel fc Hern quist 2003b, 3). 

Here, we use the 05 & G5 runs from 

ISpringel fc HernquistJ (|2003bf h which have box sizes 
of 10 & 100/i -1 Mpc, respectively. The total particle 
number is N p = 2 x 324 3 for gas and dark matter in each 
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Fig. 1. — Distribution of all DLAs (including both QSO-DLAs and GRB-host-DLAs) as a function of JVhi (panel \a]) and Ssfr (pane l 
[b]) at z = 1 — 10. The top axis in each panel indicates corresponding values of Nm and Ssfr based on the empirical Kcnnicutt (1998) 
law. 



run. The initial gas particle mass is m gas = 3.3 x 10 5 
(3.3 x 10 8 ) /i _1 Mq, and the dark matter particle mass 
is m dm = 2.1 x 1(T (2.1 x 10 9 ) hr l M Q . The comoving 
gravitational softening length, a measure of the spatial 
resolution of the simulation, is 1.2 (8.0) /i _1 kpc for 
the Q5 (G5) run. We use the Q5 run for z > 3 (for 
higher resolution), and the G5 run for z < 3 (for better 
sampling of more massive halos and longer wavelength 
perturbations). 

Previously, we have used these simulations to stud y 
the properties of DLAs (iNagamine et al.ll2004aJbl 12007), 
Ly-break galaxies at z = 3 — 6 (Nagami ne et al.l 12 004c; 
Night et al.l l2006j ) , and Lya emitters (jNagamin e et al.1 
2008). In general, the simulations show reason- 



able agreement with available galaxy observations, giv- 
ing some confidence that we are capturing the ba- 
sic aspects of hierarchical galaxy evolution in the 
context of the ACDM model. Moreover, the cos- 
mic star formation history implied by the simulations 
(jHernquist &: S pringcl 2003 ) agrees wit h observational 
estimates (jFaucher-Giguere et al.l l2008| ) and supports 
the association of GRBs wi th massive star formation 
(jFaucher-Giguere et al.1 [20081 ). The adopted cosmolog- 
ical parameters of all simulations considered here are 
(n m ,tt A ,n h ,a 8 ,h) = (0.3,0.7,0.04,0.9,0.7), where h = 
H o /(100 kms^Mpc" 1 ). 

3. RESULTS 

First, we present the distribution of all DLAs in the 
simulation as a function of Ajji in FigurefT^. The method 
of calculating iVm i n the simulations is described fully 
in iNagamine et ail (|2004al ). Briefly, we set up a cubic 
grid that covers each dark matter halo, and calculate 
Am by projecting the Hi mass distribution onto one 
of the planes. The quantity "<in" is the area cover- 
ing fraction o n the sky along the li ne element cdt (see 
Eqns. 5-7 of [Nagamin e et all [2007). and the function 
d 2 n/(d\ogN m dX) = f(N m ,X)N m ln(10) is the 'inci- 



dence rate' per unit log Ajji and per unit absorption dis- 
tance X(z), where dX — 7^y(l + z) 2 dz. The function 

f{Xm,X) is usually referred to as the column density 
distribution function. 

Figure [IK shows that, from z — 10 to z — 6, the inci- 
dence rate increases monotonically with decreasing red- 
shift at all Njji, reflecting the rapidly growing number 
of dark matter halos. From z = 6 to z = 3, there is not 
much change at log-ZVui > 22, but the number of columns 
at log JVhi < 21 has decreased significantly, which could 
owe to the UV background radiation field imposed at 
z = 6 in the simulation to model reionization. From 
z = 3 to z = 2, the number of columns at logAni > 22 
increases, but it decreases from z = 2 to z = 1, owing to 
the conversion of high-density gas into stars. 

If long GRBs are associated with the collapse of mas- 
sive stars, then their occurrence should be correlated 
with the local SFR. Therefore we define the following 
to quantify the distribution of GRB-host-DLAs: 



CgRB — host — DLA = 



d 2 n 



J SFR 



dXd log Ssfr (S S fr) 



%rb, (1) 



where Ssfr is the projected SFR in units of 
[Mq yr -1 kpc~ 2 ], and (Ssfr) is a normalization parame- 
ter. Although it is somewhat arbitrary, we take (Ssfr) 
III '. becau se it roughly cor responds to logA^i ~ 20 
based on the lKennicuttJ (|1998[ ) law. The parameter ?7grb 
denotes the mass fraction of stars that become GRBs 
and have associated afterglows for a given amount of 
star formation with a certain stellar initial mass function 
(IMF). The exact value of ?7grb depends on the IMF 
and other GRB physics. Here, we take t^grb = 10~ 3 , 
because the mass fractio n of stars with M > 8 M© 
is 23% of the total for a IChabrierl (|2003l ) IMF with a 
mass range [0.1, 100] M , and the global average of the 
GRB/SN ratio is ~ 5% rtYoon et all 12005 ISoderberd 
[20071 iCamoana et al.l 120081) . If one takes the number 
fraction (~ 6% for M > 8 M Q stars) instead, then 
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Fig. 2. — Distribution of GRB-host-DLAs as a function of 
^SFR and A r ^j HB at z = 1 — 10 is shown in the bottom panel. 
The vertical tick-marks indicate the median of the distribution for 
columns with logSsPR > -3.0 or logA^j RB > 20.3. The top 
panel show s the observed sample of 28 GRB-host-DLAs from 
IChen et al.l (120071 . broader (finer) hatching for z < 3 (z > 3)) and 
11 GRB-host-DLAs from IProchaska et all 1 120071 . yellow shade; 
excludes those at z < 1.6 with no TVhi measurements). The axis 
range of N^^ B corresponding to that of Sgpj^ is determined from 
the empirical Kcnnicutt (1998) law, divided by a factor of 2 to 
take into account of the fact that the GRB LOSs go through 
only a half of the Hi sl ab on average compared to the QSO LOSs 
(Prochaska et al. 20Q3). 

Vgrb = 3 x 10~ 4 . Here we use the mass fraction, be- 
cause weighting by £sfr/(£sfr) * s done on the basis of 
stellar mass. 

In principle, one could absorb the factor (Esfr) into 
??GRB and treat them as one parameter: ?7grb = 
7 7grb/(Ssfr), which would be the GRB rate per pro- 
jected SFR. However, here we choose to treat them sepa- 
rately to keep the physical meaning of ??grb clear. In the 
future, GRB theory may be able to estimate the value of 
?7grb ■ and observations of GRB-host-DLAs will constrain 
the ratio of t^grb / (£sfr.) (see §0]and[5]). 

It is worthwhile to look at the distribution 
d 2 n/(dXd\ogT,sFYi), before weighting it by Esfr- Fig- 
ure [1}d shows that the redshift evolution of this distribu- 
tion is stronger than in Fig. QJi. The number of columns 
with logEsFR > -2.5 (« logA^i > 21.0 for the Kenni- 
cutt law) decreases systematically from z = 1 to z = 10. 

Figure [5] shows Cgrb— host— dla as a function of 
logEsFR- Because of the weighting by £ S fr/(£sfr), 
the distribution now exhibits a peak with a broad tail at 
high Ssfr- The number of columns with logA^ni > 21 
decreases progressively from z = 1 to z = 10, owing 
to the decreasing number of massive dark matter halos 
with deep potential wells towards higher redshifts. This 
is encouraging for the use of GRB afterglows to probe 
rcionization, because the IGM attenuation signature is 
less likely to be blocked by the red damping wing of 
GRB-host-DLAs. The median value of the distribution 
at logEsFR > -3.0 (« logA^m > 20.3) is \ogN m = 
21.4, 21.1, 21.0, 20.7, 20.6, &20.4 for z = 1, 2, 3, 6, 8 &10, 



respectively. 

The distributions at z — 1 & 2 have broad peaks at 
logiVHi = 21.0 — 22.3. The current o bserved sample 
(|Chen et al.ll2007t IProchaska et al.ll2007f) is shown in the 
top panel of Figure [H and observationally there appears 
to be no indication that TVhi drops towards high-z. If 
there is any observed trend, it might even be in the op- 
posite sense, and all the observed GRB-host-DLAs are 
at z > 2. However the current observed sample is still 
small, and observational selection effects are at play. For 
example, systems with low Nm and low metal content 
are more difficult to identify, especially at higher redshift 
where the afterglows are typically fainter. We also note 
that the decline of the distribution at logEsFR < —2.5 
for z = 1 & 2 may owe to the limited resolution of the 
G5 run compared to the Q5 run. 

By integrating Cgrb— host— DLA over logS SF R, 
we obtain the 'incidence rate' of GRB-host-DLAs. 
The integral at logEsFR > —3.3 yields rates of 
(6.4, 5.7, 4.4, 1.7, 0.52, 0.14) x 10~ 4 for z = 1, 2, 3, 6, 8 &10, 
respectively, for the assumed value of ?7grb / (£sfr) = 10. 
Figure [3] compares the derived incidence rate of GRB- 
host-DLAs to that of QSO-DL As. The red data points 
are the updated version 4 of IProchaska et aTl (|2005h using 
SDSS DR5. Our simulations somewhat underpredict the 
QSO-DLA incidence rate owi ng to the underestimat e 
of f(N m ) at logjV m < 21 (jNagamine et al.l l2004al) . 
There is a stark difference between the evolution of 
the two rates: the incidence rate of GRB-host-DLAs 
decreases monotonically towards high-z, whereas the 
QSO-DLA rate increases from z = 1 to z = 6. This is 
because we assumed that GRBs are correlated with star 
formation and their distribution is not random, unlike 
the background quasars. The offset between the two 
rates tells us about the difference between the total Hi 
cross section of galaxies and the area covering fraction 
of star-forming regions. The QSO-DLA sight-lines can 
also probe the outskirts of galaxies where star formation 
is nonexistent, therefore their incidence rate is much 
higher than that of GRB-host-DLAs. 

4. PROBABILITY FOR A GIVEN GRB 

Of immediate interest to GRB observers is the chance 
probability of finding a GRB-host-DLA for a given GRB 
event. For individual GRBs, the probability of hav- 
ing a GRB-host-DLA should depend on the geometry 
of the Hi gas distribution around the GRB along the 
LOS. While our simulations do not have the resolution 
to follow the gas dynamics within molecular clouds, the 
Q5 run has a physical gravitational softening length of 
0.31 (0.18) h~ x kpc at z = 3 (6), and the gas distribu- 
tion above these scales is followed reasonably well. The 
O and B stars would have ionized all the gas within 
~ 100 pc of the GRB creating an Hn region, thus we 
expec t the DLA gas to be lo cated at > 0.1 kpc from the 
GRB (Prochas ka et al.|l2007D . Then, our simulations can 
follow the qualitative trend in the redshift evolution of 
the incidence rate of GRB-host-DLAs. 

Since our simulations r oughly match the Kennicutt law 
(|Nagamine et al.lr20 04b). it is guaranteed that columns 
with logAni > 20.3 will have star formation along the 
LOS. Therefore if long GRBs are associated with star 

4 http:/ /www. ucolick.org/~xavier/SDSSDLA/ 
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Fig. 3. — Top panel: Incidence rates of GRB-host-DLAs and 
QSO-DLAs as a function of redshift. (See text for the source 
of the data points.) For GRB-host-DLAs, the three blue curves 
correspond to »7grb/(Ssfr) = 10 2 ,10, & 1 from top to bottom. 
Bottom panel: Ratio of the two incidence rates. 



formation as we have assumed, most of the GRB LOSs 
will have a high-A^ui gas in their host galaxy, majority 
of which are GRB-host-DLAs. A more detailed analysis 
would be required to fully confirm this, by generating 
the absorption line profiles for each GRB LOS, and we 
plan to examine this in due course using higher resolution 
simulations. 

Nevertheless, Figure [5] shows that the number of high- 
Nm systems decreases with increasing redshift, and so 
we expect that the chance probability of having a high- 
JVhi DLA for a given GRB event will also decline towards 
high- z with a similar qualitative trend as shown in Fig- 
ure [3l But we stress that the incidence rate shown in 
Figure [3] is not the probability of detecting a DLA for a 
given GRB. 

5. CONCLUSIONS & DISCUSSIONS 

Using cosmological SPH simulations, we have exam- 
ined the redshift evolution of incidence rates of GRB-host- 
DLAs, assuming that long GRBs are correlated with local 
SFR. The distribution of GRB-host-DLAs is intrinsically 
different from that of QSO-DLAs, and the incidence rate 
of GRB-host-DLAs decreases monotonically towards high- 
z, whereas that of QSO-DLAs increases from z = 1 to 
z = 6. Quasars are assumed to be randomly distributed 
background sources in the sky, which illuminate the DLA 
gas in foreground galaxies. GRBs can also serve as ran- 
domly distributed beacons with respect to the DLA gas 
in foreground galaxies, but for GRB-host-DLAs, GRBs 
are not random background sources because they are in 
the same host galaxy. 



We find that the incidence rate of GRB-host-DLAs with 
logA^Hi > 21.0 decreases rapidly at z > 6, suggest- 
ing that the likelihood of observing the IGM attenua- 
tion signature in GRB afterglows increases toward higher 
rcdshifts, without being blocked by the red damping of 
DLAs in the GRB host galaxies. This enhances the 
prospects for using high-z GRBs to probe the reion- 
ization history of the Universe. Our predictions can 
be tested by upcoming high-z GRB missions, includ- 
ing JANUS (Joint Astrophysics Nascent Universe Scout) 
and SVOM (Space multi-band Variable Object Monitor). 

It might be hoped that it would be possible to esti- 
mate the incidence rate of GRB-host-DLAs by accumu- 
lating a large GRB sample. However, because GRBs arc 
not random background sources for GRB-host-DLAs, this 
would require a prohibitively large GRB sample to esti- 
mate the area covering fraction of GRB-host-DLAs from 
GRB observations alone. If long GRBs do indeed trace 
star-forming regions, and if all the LOSs to star-forming 
regions are coincident with DLAs, then one could esti- 
mate the total GRB-host-DLA incidence rate simply by 
measuring the area covering fraction of star-forming re- 
gions from deep imaging surveys of galaxies. When the 
number of GRBs becomes comparable to that of QSOs, 
one should expect non-negligible new intervening DLAs 
in GRB LOSs that are not associated with GRB hosts, 
since GRB afterglows now act as random beacons. 

An alternative possibility would be to search for 
quasars in the proximity of GRBs, or vice versa. Such 
a search of QSO-GRB pair-LOSs would yield a coin- 
cidence probability between QSO-DLAs and GRB-host- 
DLAs, which roughly corresponds to the ratio of the two 
incidence rates. For this purpose, only those systems, for 
which the QSOs are in the background and the GRBs in 
the foreground, can be used. So far, no such cases have 
been identified, but a combination of all-sky GRB sur- 
veys (e.g., BATSE, Swift, GLAST) and optical-IR imag- 
ing surveys of galaxies (e.g., SDSS, Pan-STARRS, LSST) 
may prove successful in the future. A constraint on the 
ratio of the two incidence rates would make it possible to 
estimate the ratio ?7grb/2sfr- In addition, deep obser- 
vations of GRB host galaxies could constrain Ssfr in- 
dependently. Then combining the above two constraints 
would allow us to estimate the value of t^grb • 
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